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1.  Introduction 

The  objective  of  this  study  is  to  use  a  functional  high-content  screening  (HCS)  approach  to  directly  and 
comprehensively  identify  microRNAs  (miRNAs)  that  induce  neuroblastoma  cell  differentiation  by  applying  a 
library  of  miRNA  mimics  used  to  raise  the  intracellular  levels  of  miRNAs.  Two  specific  aims  are  proposed: 
Aim  1,  identify  miRNAs  that  induce  neuroblastoma  cell  differentiation  by  HCS;  Aim  2,  Validate  the  function 
of  the  identified  miRNAs  to  induce  cell  differentiation  and  cell  growth  arrest  in  neuroblastoma  cells. 

2.  Keywords 

Neuroblastoma,  miRNA,  miRNA  mimic,  miR-449a,  neurite  outgrowth,  differentiation,  differentiation  marker, 
cell  cycle  exit,  MYCN  amplification 

3.  Overall  Project  Summary 

In  this  reporting  period,  we  mainly  focused  on  characterizing  another  novel  differentiation-inducing  miRNA, 
miR-449a.  We  demonstrated  the  differentiation-inducing  function  of  miR-449a  in  MYCN-nonamplified  and 
MYCN-amplified  neuroblastoma  cell  lines.  This  fulfills  the  objective  of  Aim  1.  We  further  validated  the  effect 
of  miR-449a  on  the  expression  of  molecular  differentiation  markers,  on  cell  cycle  distribution,  and  on  cell 
survival  and  proliferation.  This  fulfills  the  task  proposed  in  Aim  2. 

4.  Key  Research  Accomplishments 

Specific  Aim  1:  Identify  miRNAs  that  induce  neuroblastoma  cell  differentiation  by  HCS. 

The  objective  of  this  Aim  is  to  identify  miRNA  mimics  that  potentially  have  general  differentiation- inducing 
effect  in  neuroblastoma  cell  lines  regardless  of  the  genetic  backgrounds  of  the  cell  lines.  We  have  completed  the 
screen  in  the  MYCN-amplified  cell  line  BE(2)-C  cells,  and  we  have  identified  a  group  of  novel  miRNA  mimics 
that  induce  BE(2)-C  cell  differentiation.  Although  we  have  not  completed  the  screen  in  the  MYCN- 
nonamplified  cells,  we  did  demonstrate  that  several  differentiation-inducing  miRNAs  identified  in  BE(2)-C 
cells  also  induce  differentiation  in  MYCN-nonamplified  cell  lines,  and  we  have  reported  these  miRNAs  in  the 
Progress  Report  in  the  previous  grant  period. 

Here,  we  report  an  additional  miRNA,  miR-449a  that  have  general  differentiation-inducing  effect  in 
neuroblastoma  cell  lines  with  distinct  genetic  backgrounds,  including  the  MYCN  amplification  status. 

4.1.  miR-449a  overexpression  induces  cells  differentiation  in  neuroblastoma  cells  with  different  genetic 
backgrounds. 

In  our  HCS,  we  identified  that  miR-449a  mimic  significantly  induces  neurite  outgrowth  in  BE(2)-C  cells  at 
25nM  (1).  To  further  characterize  the  differentiation- inducing  function  of  miR-449a,  we  examined  the  time- 
dependent  and  dose-dependent  effect  of  miR-449a  mimic  on  neurite  outgrowth  in  BE(2)-C  cells.  As  shown  in 
Figure  1A-B,  miR-449a  mimic  induces  neurite  elongation  in  both  time-  and  dose-dependent  manners,  and  the 
effect  is  significant  at  very  low  concentrations  of  miR-449a  mimic  (as  low  as  0.5  nM).  We  further  confirmed 
the  effect  of  miR-449a  on  neurite  outgrowth  using  a  miR-449a  precursor  mimic.  A  miRNA  precursor  mimic  is  a 
partially  double-stranded  RNA  designed  to  mimic  the  endogenous  miRNA  precursor,  in  contrast  to  the  fully 
complementary  double-stranded  design  of  a  miRNA  mimic  that  recapitulates  the  mature  miRNA.  Figure  1C-D 
shows  that  the  effect  of  the  miR-449a  precursor  mimic  on  neurite  outgrowth  is  comparable  to  that  of  the  miR- 
449a  mimic,  indicating  that  the  differentiation-inducing  function  of  miR-449a  mimic  is  unlikely  to  be  caused  by 
off-target  effects. 

In  order  to  examine  whether  the  differentiation-inducing  effect  of  miR-449a  mimic  in  neuroblastoma  is 
universal,  we  further  investigated  the  effect  of  miR-449a  mimic  on  cell  differentiation  in  multiple 
neuroblastoma  cell  lines  with  different  genetic  backgrounds  (Table  1).  As  shown  in  Figure  1E-F,  miR-449a 
mimic  significantly  induces  neurite  outgrowth  in  all  of  these  neuroblastoma  cells  as  compared  to  control, 
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although  the  potency  of  its  effect  distinctly  varies  among  cell  lines.  These  results  support  the  general 
differentiation-inducing  function  of  miR-449a  in  neuroblastoma  cells. 


f  miR-449a  precursor 

Control  oligo  miR-449a  mimic  mimic 


oligo  mimic  precursor 


Figure  1.  miR-449a  functions  as  an  inducer  of  cell  differentiation  in  neuroblastoma  cells.  (A-B),  Time-  (A)  and  Dose-  (B)  dependent 
effect  of  miR-449a  mimic  on  neurite  outgrowth  in  BE(2)-C  cells.  2500  cells  were  transfected  with  different  concentrations  of  miR-449a 
mimic  or  negative  control  oligo  in  triplicates  in  96-well  plates,  and  neurite  outgrowth  was  measured  every  6  hours  for  4  days.  Relative 
neurite  length  is  defined  as  neurite  length  per  cell  body  area.  Shown  are  the  normalized  relative  neurite  length,  presented  as  average  values 
of  the  three  replicates.  *,/?<0.05,  comparing  to  the  negative  control  oligo  at  the  corresponding  time  or  dose  points.  (C-D),  Effect  of  miR- 
449a  mimic  and  precursor  mimic  on  neurite  outgrowth  in  BE(2)-C  cells.  Cells  were  transfected  with  25  nM  control  oligo,  miR-449a  mimic 
or  precursor  mimic  in  triplicates,  with  neurite  lengths  measured  as  above  after  4  days.  Shown  are  the  representative  cell  images  (C) 
analyzed  to  define  neurites  (pink)  and  cell  body  areas  (yellow),  and  the  quantification  of  neurite  length  (D).  *,  /?<0.05,  comparing  to 
control.  (E-F),  Effect  of  miR-449a  mimic  on  neurite  outgrowth  in  multiple  neuroblastoma  cell  lines.  Cells  were  transfected  with  25  nM 
miR-449a  mimic  or  control  oligo  in  triplicates.  Neurite  lengths  were  measured  as  above.  E,  Representative  cell  images  analyzed  to  define 
neurites  (pink)  and  cell  body  areas  (yellow)  after  4  days  of  transfection.  F,  Quantification  of  neurite  length.  *,  /?<0.05,  comparing  to 
control. 


Cell  line 

Age 

Gender 

Stage 

Chr  lp  alteration 

Chr  17  alteration 

MYCN  ampl. 

BE(2)-C 

2 

male 

4 

Y 

N 

Y 

SK-N-BE 

1.8 

male 

4 

Y 

N 

Y 

LAN6 

5.8 

male 

4 

N 

unk 

N 

BE(2)-M17 

2 

male 

4 

Y 

N 

Y 

KELLY 

unk 

unk 

unk 

unk 

unk 

Y 

Table  1.  Genetic  backgrounds  of 
neuroblastoma  cell  lines  used  in  the 
study.  Shown  are  the  name  of  the  cell 
line,  age  and  gender  of  the  patient, 
stage  of  the  tumor  from  which  the  cell 
lines  are  derived,  chromosome  lp  and 
17  alterations,  and  MYCN  gene 
amplification  status.  *unk,  unknown; 
Chr,  Chromosome;  ampl,  amplification. 
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Specific  Aim  2.  Validate  the  function  of  the  identified  miRNAs  in  inducing  cell  differentiation  and  cell 
growth  arrest  in  neuroblastoma  cell  lines. 

The  goal  of  this  Aim  is  to  use  2nd  and  3rd  tier  analyses  to  validate  the  candidate  miRNA  mimics  identified 
from  the  HCS,  by  examining  the  effects  of  the  miRNA  mimics  on  inducing  neuroblastoma  cell  differentiation 
and  cell  growth  arrest.  Focusing  on  miR-449a,  we  validated  its  tumor  suppressive  function  by  examining  its 
effect  on  expression  of  molecular  differentiation  markers,  on  cell  survival  and  proliferation,  and  on  cell  cycle 
distribution. 

Aim  2 A.  Validate  the  function  of  the  identified  miRNAs  in  inducing  cell  differentiation  by  examining  the 
expression  of  molecular  differentiation  markers. 

4.2.  miR-449a  mimic  induces  expression  of  molecular  differentiation  markers  in  multiple  neuroblastoma 
cell  lines. 

In  order  to  further  validate  the  differentiation-inducing  function  of  miR-449a,  we  examined  the  effect  of 
miR-449a  mimic  on  the  expression  of  molecular  differentiation  markers  in  multiple  neuroblastoma  cell  lines. 
Consistent  with  its  effect  on  neurite  outgrowth,  as  shown  in  Figure  2,  miR-449a  mimic  significantly  increases 
the  expression  of  neuroblastoma  cell  differentiation  markers,  including  pill-tubulin,  neuron  specific  enolase 
(NSE)  and  growth  associated  protein  43  (GAP43)  (2-4),  in  all  the  five  cell  lines  tested. 
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Figure  2.  Effect  of  miR-449a  mimic  on  protein  expression  levels  of  differentiation  markers  in  multiple  neuroblastoma  cell  lines.  Cells  were 
treated  with  25  nM  miR-449a  mimic  or  control  oligos  (Control)  for  4  days.  Expression  of  differentiation  markers  pill-tubulin,  GAP43  and  NSE 
were  determined  by  Western  blots  with  calnexin  protein  levels  used  as  a  loading  control.  Shown  are  (A)  the  images  of  the  Western  blots  as  well  as 
(B-F)  the  quantification  of  the  band  intensities  using  Image  J.  The  band  intensities  associated  with  the  differentiation  markers  were  normalized  to 
those  associated  with  calnexin  to  allow  comparison  between  samples.  Shown  are  the  average  results  from  three  independent  experiments.  *, 
p<0.05,  comparing  to  Control. 
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Aim  2b.  Validate  the  function  of  the  identified  miRNA  mimics  in  inducing  cell  growth  arrest  of 
neuroblastoma  cells  by  cell  cycle  analysis. 

4.3.  Overexpression  of  miR-449a  induces  G0/G1  arrest  and  cell  cycle  exit  in  neuroblastoma  cells. 
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Figure  3.  Effect  of  miR-449a  on  cell  cycle  progression  in  neuroblastoma  cells.  A-C,  Effect  of  miR-449a  overexpression  on  cell  cycle 
distribution  in  (A)  BE(2)-C,  (B)  SKNBE  and  (C)  KELLY  cells.  Cells  were  transfected  with  25  nM  miR-449a  mimic  or  control  oligo,  or 
treated  with  5  pM  RA.  Mock  transfected  cells  were  used  as  an  additional  negative  control.  After  3  days,  cells  were  collected  and  stained  with 
propidium  iodide  for  cell  cycle  analysis.  The  fraction  of  cells  in  G0/G1,  S  and  G2  phases  was  quantified  using  the  Jean-Jett-Fox  model. 
Similar  results  were  obtained  from  three  independent  experiments.  D-G,  Effect  of  miR-449a  and  RA  on  DNA  synthesis  as  measured  by 
BrDU  incorporation  in  BE(2)-C  (D),  SKNBE  (E)  and  KELLY  (F)  cells.  Cells  were  transfected  with  miR-449a  mimic  or  control  oligos  (25 
nM),  or  treated  with  RA  (5  pM)  for  4  days.  Cells  were  incubated  with  BrdU  during  the  final  2  h  of  culture.  Cells  were  then  stained  for  BrDU 
incorporation  (red)  as  described  in  the  Materials  and  Methods.  Cell  nuclei  were  identified  by  nuclear  staining  with  DAPI  (blue).  G, 
Quantification  of  percentage  of  cells  positive  for  BrDU  staining.  *,/?<0.05,  comparing  to  control. 
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Cell  cycle  arrest  is  another  marker  that  indicates  the  differentiation  of  neuroblastoma  cells  (5).  We  therefore 
further  examined  the  effect  of  miR-449a  mimic  on  cell  cycle  distribution  and  on  cell  cycle  exit.  As  shown  in 
Figure  3A-C,  comparing  to  the  control  cells,  miR-449a  overexpression  leads  to  marked  increase  of  cell 
numbers  in  G0/G1  phase  and  dramatic  decrease  of  cells  numbers  in  S-  and  G2/M  phases,  suggesting  that  miR- 
449a  blocks  cell  cycle  progression  through  arresting  cells  in  G0/G1  phase.  We  further  show  that  the  effect  of 
miR-449a  on  cell  cycle  distribution  is  comparable  to  that  induced  by  RA,  which  is  a  differentiation  agent  that  is 
currently  used  as  standard  of  care  for  high-neuroblastoma  therapy.  It  has  been  known  that  RA  treatment  of 
neuroblastoma  cells  lead  to  terminal  cell  differentiation,  which  is  accompanied  with  cell  cycle  exit  (6).  We 
therefore  examined  whether  miR-449a  overexpression  has  similar  effect  in  neuroblastoma  cells  by  examining 
the  BrDU  incorporation  into  DNA  in  the  nucleus  of  individual  cells  using  immunofluorescence  staining.  As 
shown  in  Figure  3D-G,  comparing  to  the  control  cells,  more  cells  treated  with  miR-449a  mimic  or  RA  are 
negative  for  BrDU  staining,  which  indicates  that  these  cells  failed  to  enter  S-phase.  These  results  indicate  that 
induction  of  cell  cycle  exit  is  achieved  by  miR-449a  overexpression  in  neuroblastoma  cells,  and  the  effect  of 
miR-449a  in  inducing  cell  cycle  exit  is  comparable  to  RA. 

4.4.  miR-449a  reduces  neuroblastoma  cell  survival  and  proliferation. 

To  further  characterize  the  potential  tumor  suppressive  function  of  miR-449a  in  neuroblastoma,  we  examined 
whether  the  induced  cell  differentiation  by  miR-449a  is  coupled  with  reduced  neuroblastoma  cell  survival.  As 
show  in  Figure  4A,  miR-449a  mimic  decreases  cell  viability  in  a  dose-dependent  manner  in  all  the  examined 
neuroblastoma  cell  lines.  We  further  investigated  the  mechanisms  underlying  the  reduced  cell  survival  induced 
by  miR-449a.  Figure  4B  shows  that  miR-449a  mimic  activates  the  apoptotic  pathway  in  neuroblastoma  cells,  as 
measured  by  increased  levels  of  cleaved  caspase  3.  To  further  demonstrate  its  tumor  suppressive  function,  we 
next  examined  the  impact  of  the  miR-449a  mimic  on  the  tumorigenic  potential  of  neuroblastoma  cells  by  colony 
formation  assay.  As  shown  in  Figure  4C-D,  miR-449a  mimic  significantly  decreases  the  capacity  of  the 
neuroblastoma  cells  to  form  colonies,  indicating  that  miR-449a  reduces  the  tumorigenicity  of  neuroblastoma 
cells.  Altogether,  our  results  support  that  miR-449a  functions  as  a  strong  suppressor  of  tumorigenicity  in 
neuroblastoma,  and  this  is  achieved  by  inducing  cell  differentiation,  reducing  cell  proliferation  and  activating 
the  cell  apoptotic  pathway. 

/\  g  be(2)-c  kelly  Figure  4.  Effect  of  miR-449a  mimic  on 

neuroblastoma  cell  survival  and  proliferation.  A, 

Dose-dependent  effect  of  miR-449a  mimic  on  cell 
viability  in  multiple  neuroblastoma  cell  lines.  Cells 
were  transfected  with  different  concentrations  of 
miR-449a  mimic  in  triplicates  in  96-well  plates. 
After  4  days,  cell  viability  was  measured  using  the 
CellTiter-Glo  Luminescent  Cell  Viability  Assay. 
Shown  are  the  relative  cell  viabilities,  presented  as 
averages  of  the  three  replicates,  normalized  to  the 
viabilities  of  cells  at  the  0  nM  dose  point  (mock 
transfected  cells).  *,/?<0.05,  comparing  to  the  0  nM 
dose  point.  B,  Western  blot  analysis  of  Caspase-3 
activation.  Cells  were  transfected  with  25  nM  miR- 
449a  mimic  or  control  oligos  for  72  h,  after  which 
cell  lysates  were  harvested.  Intact  and  cleave 
caspase-3  was  detected  using  an  anti-caspase-3 
antibody  that  recognizes  both  the  intact  and  the 
cleaved  caspase-3.  Calnexin  levels  were  measured 
as  a  loading  control.  (C-D)  Colony  formation  assay 
as  a  function  of  miR-449a  mimic.  Treatments  were 
conducted  in  five  replicates.  Show  are  the 
representative  images  (C)  and  quantification  of  the 
number  of  colonies  (D).  *  /K0.05,  comparing  to 
Control. 
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5.  Conclusion  and  future  plans 

In  summary,  we  have  demonstrated  that  miR-449a  functions  as  an  inducer  of  cell  differentiation  in 
neuroblastoma  cell  lines  with  distinct  genetic  backgrounds,  including  the  MYCN  amplification  status.  We 
further  show  that  miR-449s  mimic  induces  cell  cycle  arrest  and  cell  cycle  exit,  induces  cell  apoptosis,  and 
decreases  the  capacity  of  the  neuroblastoma  cells  to  form  colonies.  All  these  results  support  the  strong  tumor 
suppressive  function  of  miR-449a  in  neuroblastoma  cells,  providing  the  experimental  evidence  for  developing 
miR-499a-based  differentiation  therapy  for  neuroblastoma. 

In  the  next  grant  period,  we  will  (1)  complete  the  HCS  screen  in  MYCN-nonamplified  cells  and  identify 
additional  miRNAs  that  potentially  have  general  effect  on  cell  differentiation  in  neuroblastoma,  and  (2)  validate 
the  rest  of  the  identified  differentiation-inducing  miRNAs  by  examining  expression  of  the  molecular 
differentiation  markers,  cell  survival  and  cell  cycle  distribution. 

6.  Publications,  Abstracts  and  Presentations 

6.1.  Publications 

Zhao  Z,  Ma  X,  Sung  D,  Li  M,  Kosti  A,  Lin  G,  Chen  Y,  Pertsemlidis  A,  Hsiao  TH  and  *Du  L.  microRNA- 
449a  functions  as  a  tumor  suppressor  in  neuroblastoma  through  inducing  cell  differentiation  and  cell  cycle 
arrest.  RNA  Biology.  2015;  12(5):  538-554.  *,  Corresponding  author. 

6.2.  Abstracts 

1.  Ma  X,  Li  M,  Zhao  Z,  Pertsemlidis  A,  Sung  D  and  Du  L.  Crosstalk  between  MYCN  and  differentiation- 
inducing  microRNAs  in  neuroblastoma.  Greehey  Children ’s  Cancer  Research  Institute  2014  Symposium. 
University  of  Texas  Health  Science  Center  at  San  Antonio.  San  Antonio,  TX.  November  21,  2014. 

2.  Zhao  Z,  Ma  X,  Li  M,  Kosti  A,  Lin  G,  Chen  Y,  Pertsemlidis  A,  Hsiao  TH  and  Du  L.  microRNA-449a 
functions  as  a  tumor  suppressor  in  neuroblastoma  through  inducing  cell  differentiation  and  cell  cycle  arrest. 
Greehey  Children’s  Cancer  Research  Institute  2014  Symposium.  University  of  Texas  Health  Science  Center 
at  San  Antonio.  San  Antonio,  TX.  November  21,  2014. 

6.3.  Presentations 

1.  Identifying  novel  differentiation  agents  for  neuroblastoma  therapy.  Innovative  Drug  Discovery  & 
Nanotechnology  session.  2015  Drug  Discovery  and  Therapy  World  Congress.  Boston,  USA.  July  22- 
25th,  2015. 

2.  Using  high-content  screening  to  identify  novel  differentiation  agents  for  neuroblastoma.  Session  603: 
Lead  Discovery  and  Optimization.  2015  BIT‘s  8th  Annual  World  Cancer  Congress.  Beijing,  China.  May 
15-17th,  2015. 

3.  Targeting  the  differentiation  pathway  in  neuroblastoma  differentiation  therapy.  Department  of 
Chemistry  and  Biochemistry.  Texas  State  University.  San  Marcos,  TX.  December  04,  2014. 

4.  Differentiation  therapy  in  neuroblastoma.  2014  EDT  Program  Retreat.  Cancer  Therapy  and  Research 
Center,  University  of  Texas  Health  Science  Center  at  San  Antonio.  San  Antonio,  TX.  November  7, 
2014. 

7.  Inventions,  Patents  and  Licenses 

None  in  this  grant  period. 

8.  Reportable  Outcome 

8.1.  We  demonstrated  that  our  identified  novel  differentiation- inducing  miRNA-449a  has  general  effect  in 
inducing  cell  differentiation  in  multiple  neuroblastoma  cell  lines  with  distinct  genetic  backgrounds,  including 
MYCN-nonamplified  and  MYCN-amplified  neuroblastoma  cells. 

8.2.  We  demonstrated  the  miR-449a  mimic  induces  neuroblastoma  cell  cycle  exit,  reduces  cell  viability  and 
reduces  the  capacity  of  neuroblastoma  cells  to  form  colonies,  suggesting  that  miR-449a  suppresses 
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tumorigenecity  in  neuroblastoma.  This  evidence  strongly  supports  the  therapeutic  potential  of  miR-449a  mimic 
in  neuroblastoma. 

9.  Other  Achievements 

Based  on  the  work  supported  by  this  award,  I  submitted  the  following  grant  applications: 

1)  NIH/NCI  (1R01CA194182-01A1)  (R01)  (L.  Du)  12/01/2015-11/30/2020  4.8  Cal.  mos. 

($1250,000  direct/$1853,900  total) 

Therapeutic  mechanisms  of  differentiation-inducing  microRNAs  in  neuroblastoma 

Goal:  To  define  the  role  of  microRNAs  in  regulating  neuroblastoma  cell  differentiation,  and  to  develop 

miRNA-based  differentiation  agents  for  treating  neuroblastoma. 

Role:  Principal  Investigator 
Status:  pending 

2)  Pew  Charitable  Trust  (L.  Du)  07/01/2015-06/30/2019  2.4  Cal.  mos. 

($240,000  direct) 

Differentiation-inducing  microRNAs  in  neuroblastoma  therapeutics 

Goal:  To  develop  novel  microRNA-based  differentiation  therapies  for  treating  neuroblastoma. 

Role:  Principal  Investigator 
Status:  pending 
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microRNA-449a  (miR-449a)  has  been  identified  to  function  as  a  tumor  suppressor  in  several  types  of  cancers. 
However,  the  role  of  miR-449a  in  neuroblastoma  has  not  been  intensively  investigated.  We  recently  found  that  the 
overexpression  of  miR-449a  significantly  induces  neuroblastoma  cell  differentiation,  suggesting  its  potential  tumor 
suppressor  function  in  neuroblastoma.  In  this  study,  we  further  investigated  the  mechanisms  underlying  the  tumor 
suppressive  function  of  miR-449a  in  neuroblastoma.  We  observed  that  miR-449a  inhibits  neuroblastoma  cell  survival 
and  growth  through  2  mechanisms — inducing  cell  differentiation  and  cell  cycle  arrest.  Our  comprehensive 
investigations  on  the  dissection  of  the  target  genes  of  miR-449a  revealed  that  3  novel  targets-  MFAP4,  PKP4  and 
TSEN15  -play  important  roles  in  mediating  its  differentiation-inducing  function.  In  addition,  we  further  found  that  its 
function  in  inducing  cell  cycle  arrest  involves  down-regulating  its  direct  targets  CDK6  and  LEF1.  To  determine  the 
clinical  significance  of  the  miR-449a-mediated  tumor  suppressive  mechanism,  we  examined  the  correlation  between 
the  expression  of  these  5  target  genes  in  neuroblastoma  tumor  specimens  and  the  survival  of  neuroblastoma  patients. 
Remarkably,  we  noted  that  high  tumor  expression  levels  of  all  the  3  miR-449a  target  genes  involved  in  regulating  cell 
differentiation,  but  not  the  target  genes  involved  in  regulating  cell  cycle,  are  significantly  correlated  with  poor  survival 
of  neuroblastoma  patients.  These  results  suggest  the  critical  role  of  the  differentiation-inducing  function  of  miR-449a  in 
determining  neuroblastoma  progression.  Overall,  our  study  provides  the  first  comprehensive  characterization  of  the 
tumor-suppressive  function  of  miR-449a  in  neuroblastoma,  and  reveals  the  potential  clinical  significance  of  the  miR- 
449a-mediated  tumor  suppressive  pathway  in  neuroblastoma  prognosis. 


Introduction 

Neuroblastoma  is  one  of  the  most  common  solid  cancers  of 
childhood.  High-risk  neuroblastoma  is  one  of  the  leading  causes 
of  cancer-related  deaths  in  childhood,1,2  because  only  a  few  high- 
risk  neuroblastoma  patients  become  long-term  survivors  with  cur¬ 
rently  available  therapeutic  agents  for  treating  neuroblastoma.  Dif¬ 
ferentiation  therapy  was  developed  based  on  the  knowledge  that 
neuroblastoma  arises  from  the  neural  crest  cell  precursors  that  fail 
to  complete  the  differentiation  process.2,3  It  is  an  approach  to 


induce  malignant  cells  to  differentiate  into  mature  cells  and 
thereby  leading  to  tumor  growth  arrest.2,4-6  Currently,  the  differ¬ 
entiation  agent  1 3-cis- retinoic  acid  (RA)  is  used  as  the  standard  of 
care  for  post-remission  maintenance  therapy  in  high-risk  neuro¬ 
blastoma,2  highlighting  the  critical  role  of  differentiation  therapy 
in  neuroblastoma  treatment.  However,  resistance  to  RA  treatment 
is  common,  and  more  than  50%  of  patients  treated  with  RA  still 
subsequently  develop  recurrent  neuroblastoma.7,8  The  develop¬ 
ment  of  new  differentiation  agents  has  been  greatly  hampered  due 
to  the  poor  understanding  of  the  mechanisms  that  control 
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neuroblastoma  cell  differentiation.  This  highlights  the  need  to  fur¬ 
ther  investigate  the  mechanisms  that  regulate  neuroblastoma  cell 
differentiation  in  order  to  identify  novel  therapeutic  targets  for 
developing  new  and  more  effective  differentiation  therapies. 

microRNAs  (miRNAs)  are  endogenously  expressed  small 
RNAs  that  play  a  critical  role  in  tumorigenesis,9"13  and  miRNA- 
based  anti-cancer  therapeutics  are  under  intensive  investiga¬ 
tion.14"17  One  of  the  most  promising  approaches  that  have  been 
developed  is  called  miRNA  replacement  therapy,  in  which  the 
intracellular  levels  of  tumor  suppressive  miRNAs  are  raised  using 
miRNA  mimics,  synthetic  oligonucleotides  that  mimic  the 
sequences  and  function  of  the  endogenous  miRNAs.18  These 
miRNA  mimics  have  been  used  to  effectively  treat  cancers  in  ani¬ 
mal  studies,14’15  showing  the  promise  of  developing  them  into 
anti-cancer  agents.14"17  Recent  studies  have  shown  that  miRNAs 
play  an  important  role  in  regulating  neuroblastoma  cell  differen¬ 
tiation,19"25  thus  suggesting  the  potential  application  of  miRNA- 
based  therapeutics  in  neuroblastoma  differentiation  therapy. 

To  this  end,  we  recently  identified  a  group  of  novel  miRNAs 
that  function  as  potent  inducers  of  neuroblastoma  cell  differenti¬ 
ation,20  which  provides  candidates  for  developing  novel  miRNA 
replacement  therapies  in  neuroblastoma.  One  of  the  most  potent 
differentiation-inducing  miRNAs  that  we  identified  was  miR- 
449a,  which  has  been  demonstrated  to  function  as  an  important 
tumor  suppressor  in  several  types  of  cancers.26"33  This  evidence 
highlights  the  universal  tumor  suppressive  function  of  miR-449a 
in  cancers  and  supports  the  potential  application  of  miR-449a 
replacement  therapy  in  cancer  treatment.  In  neuroblastoma, 
miR-449a  has  been  shown  to  target  the  3'  untranslated  region 
(3rUTR)  of  the  MYCN  mRNA.34  However,  the  tumor  suppres¬ 
sive  function  of  miR-449a  in  neuroblastoma  has  not  been  fully 
evaluated,  and  the  molecular  mechanisms  underlying  its  tumor 
suppressive  function  as  well  as  the  clinical  relevance  of  this  tumor 
suppressive  mechanism  is  not  defined.  In  the  current  study,  we 
aim  to  further  define  the  tumor  suppressive  function  of  miR- 
449a  in  neuroblastoma,  by  elucidating  the  underlying  mecha¬ 
nisms  for  such  function  and  evaluating  the  potential  clinical  sig¬ 
nificance  of  the  miR-449a-mediated  tumor  suppressive  pathway 
in  neuroblastoma  development. 

Results 

miR-449a  overexpression  induces  cells  differentiation  in 
neuroblastoma  cells  with  different  genetic  backgrounds 

We  previously  found  that  increasing  intracellular  levels  of 
miR-449a  with  a  miR-449a  mimic  significantly  induces  cell  dif¬ 
ferentiation  in  a  neuroblastoma  cell  line  BE(2)-C.20  Here,  to  fur¬ 
ther  characterize  the  differentiation-inducing  function  of  miR- 
449a  in  neuroblastoma,  we  first  examined  the  time-dependent 
and  dose-dependent  effect  of  miR-449a  mimic  on  neurite  out¬ 
growth  in  BE(2)-C  cells.  As  shown  in  Figure  1A-B,  miR-449a 
mimic  induces  neurite  elongation  in  both  time-  and  dose-depen¬ 
dent  manners,  and  the  effect  is  significant  at  very  low  concentra¬ 
tions  of  miR-449a  mimic  (as  low  as  0.5  nM).  We  further 
confirmed  the  effect  of  miR-449a  on  neurite  outgrowth  using  a 


miR-449a  precursor  mimic.  A  miRNA  precursor  mimic  is  a  par¬ 
tially  double-stranded  RNA  designed  to  mimic  the  endogenous 
miRNA  precursor,  in  contrast  to  the  fully  complementary  dou¬ 
ble-stranded  design  of  a  miRNA  mimic  that  recapitulates  the 
mature  miRNA.  Figure  1C-D  shows  that  the  effect  of  the  miR- 
449a  precursor  mimic  on  neurite  outgrowth  is  comparable  to 
that  of  the  miR-449a  mimic,  indicating  that  the  differentiation- 
inducing  function  of  miR-449a  mimic  is  unlikely  to  be  caused 
by  off-target  effects.  We  further  investigated  the  effect  of  miR- 
449a  mimic  on  cell  differentiation  in  multiple  neuroblastoma 
cell  lines  with  different  genetic  backgrounds  (Table  SI).  As 
shown  in  Figure  IE— F,  miR-449a  mimic  significantly  induces 
neurite  outgrowth  in  all  of  these  neuroblastoma  cells  as  compared 
to  control,  although  the  potency  of  its  effect  distinctly  varies 
among  cell  lines.  Consistent  with  its  effect  on  neurite  outgrowth, 
as  shown  in  Figure  1G  and  Supplementary  Figure  1,  miR-449a 
mimic  significantly  increases  the  expression  of  neuroblastoma 
cell  differentiation  markers,  including  |3III-tubulin,  neuron  spe¬ 
cific  enolase  (NSE)  and  growth  associated  protein  43 
(GAP43), 35-37  in  all  the  5  cell  lines  tested. 

Next  we  were  interested  in  examining  whether  the  endoge¬ 
nous  expression  of  miR-449a  is  regulated  during  neuroblastoma 
cell  differentiation.  We  measured  the  endogenous  miR-449a  lev¬ 
els  in  BE(2)-C  cells  that  were  induced  into  differentiation  by  RA 
(Fig.  S2).  As  shown  in  Figure  1H,  the  endogenous  expression  of 
miR-449a  in  differentiated  BE(2)-C  cells  is  significantly 
increased  comparing  to  the  undifferentiated  (Control)  cells. 
These  results  strongly  suggest  that  the  endogenous  expression  of 
miR-449a  is  suppressed  in  undifferentiated  neuroblastoma  cells. 

miR-449a  overexpression  reduces  cell  proliferation  and 
induces  apoptosis  in  neuroblastoma  cells 

To  further  characterize  the  potential  tumor  suppressive  func¬ 
tion  of  miR-449a  in  neuroblastoma,  we  examined  whether  the 
induced  cell  differentiation  by  miR-449a  is  coupled  with  reduced 
neuroblastoma  cell  survival.  As  shown  in  Figure  2A,  miR-449a 
mimic  decreases  cell  viability  in  a  dose-dependent  manner  in  all 
the  examined  neuroblastoma  cell  lines.  We  further  investigated 
the  mechanisms  underlying  the  reduced  cell  survival  induced  by 
miR-449a.  As  shown  in  Figure  2B,  miR-449a  overexpression  sig¬ 
nificantly  decreases  DNA  synthesis  in  the  cells,  as  measured  by  the 
reduced  bromodeoxyuridine  (BrdU)  incorporation  into  cell 
DNA,  indicating  its  direct  effect  in  inhibiting  cell  proliferation. 
Figure  2C  shows  that  miR-449a  mimic  also  activates  the  apopto- 
tic  pathway  in  neuroblastoma  cells,  as  measured  by  increased  levels 
of  cleaved  caspase  3.  To  further  demonstrate  its  tumor  suppressive 
function,  we  next  examined  the  impact  of  the  miR-449a  mimic 
on  the  tumorigenic  potential  of  neuroblastoma  cells  by  colony  for¬ 
mation  assay.  As  shown  in  Figure  2D— E,  miR-449a  mimic  signif¬ 
icantly  decreases  the  capacity  of  the  neuroblastoma  cells  to  form 
colonies,  indicating  that  miR-449a  reduces  the  tumorigenicity  of 
neuroblastoma  cells. 

Altogether,  our  results  support  that  miR-449a  functions  as  a 
strong  suppressor  of  tumorigenicity  in  neuroblastoma,  and  this  is 
achieved  by  inducing  cell  differentiation,  reducing  cell  prolifera¬ 
tion  and  activating  the  cell  apoptotic  pathway. 
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Figure  1.  miR-449a  functions  as  an  inducer  of  cell  differentiation  in  neuroblastoma  cells.  (A-B),  Time-  (A)  and  Dose-  (B)  dependent  effect  of  miR- 
449a  mimic  on  neurite  outgrowth  in  BE(2)-C  cells.  2500  cells  were  transfected  with  different  concentrations  of  miR-449a  mimic  or  negative  control  oligo 
in  triplicates  in  96-well  plates,  and  neurite  outgrowth  was  measured  every  6  h  for  4  d.  Relative  neurite  length  is  defined  as  neurite  length  per  cell  body 
area.  Shown  are  the  normalized  relative  neurite  lengths,  presented  as  average  values  of  the  3  replicates.  *,  P  <  0.05,  comparing  to  the  negative  control 
oligo  at  the  corresponding  time  or  dose  points.  (C-D),  Effect  of  miR-449a  mimic  and  precursor  mimic  on  neurite  outgrowth  in  BE(2)-C  cells.  Cells  were 
transfected  with  25  nM  control  oligo,  miR-449a  mimic  or  precursor  mimic  in  triplicates,  with  neurite  lengths  measured  as  above  after  4  d.  Shown  are  the 
representative  cell  images  (C)  analyzed  to  define  neurites  (pink)  and  cell  body  areas  (yellow),  and  the  quantification  of  neurite  length  (D).  *,  P  <  0.05, 
comparing  to  control.  (E-F),  Effect  of  miR-449a  mimic  on  neurite  outgrowth  in  multiple  neuroblastoma  cell  lines.  Cells  were  transfected  with  25  nM  miR- 
449a  mimic  or  control  oligo  in  triplicates.  Neurite  lengths  were  measured  as  above.  (E)  Representative  cell  images  analyzed  to  define  neurites  (pink)  and 
cell  body  areas  (yellow)  after  4  d  of  transfection.  (F)  Quantification  of  neurite  lengths.  *,  P  <  0.05,  comparing  to  control.  G,  Effect  of  miR-449a  overexpres¬ 
sion  on  the  protein  expression  levels  of  cell  differentiation  markers  (3111-tubulin,  NSE  and  GAP43  with  calnexin  protein  levels  used  as  a  loading  control. 
Cells  were  transfected  with  25  nM  miR-449a  mimic  or  control  oligo,  and  protein  levels  were  determined  by  Western  blot  after  4  d.  (H)  Endogenous 
expression  levels  of  miR-449a  in  undifferentiated  and  differentiated  BE(2)-C  cells.  BE(2)-C  cells  were  treated  with  10  |xM  of  RA  or  the  carrier  DMSO  (Con¬ 
trol)  for  5  d  to  induce  cell  differentiation.  RNA  was  then  isolated,  and  expression  of  miR-449a  in  cells  were  measured  by  qPCR  with  levels  of  1 8s  rRNA  as  a 
loading  control.  *,  P  <  0.05,  comparing  to  control. 
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Figure  2.  Effect  of  miR-449a  mimic  on  neuroblastoma  cell  survival  and  proliferation.  (A)  Dose-dependent  effect  of 
miR-449a  mimic  on  cell  viability  in  multiple  neuroblastoma  cell  lines.  Cells  were  transfected  with  different  concentrations 
of  miR-449a  mimic  in  triplicates  in  96-well  plates.  After  4  d,  cell  viability  was  measured  using  the  CellTiter-Glo  Luminescent 
Cell  Viability  Assay.  Shown  are  the  relative  cell  viabilities,  presented  as  averages  of  the  3  replicates,  normalized  to  the  via¬ 
bilities  of  cells  at  the  0  nM  dose  point  (mock  transfected  cells).  *,  P  <  0.05,  comparing  to  the  0  nM  dose  point.  (B)  BrdU 
incorporation  as  a  function  of  miR-449a  overexpression.  Cells  were  transfected  with  25  nM  miR-449a  mimic  or  control  oli- 
gos  (Control).  After  4  d,  BrdU  incorporation  was  determined  as  described  in  the  Materials  and  Methods.  Shown  are  the  rel¬ 
ative  BrDU  incorporation  levels  with  the  BrDU  levels  associated  with  cells  treated  with  miR-449a  normalized  to  those 
treated  with  control  oligos  for  each  of  the  cell  lines.  *  P  <  0.05,  comparing  to  control.  (C)  Western  blot  analysis  of  Cas- 
pase-3  activation.  Cells  were  transfected  with  25  nM  miR-449a  mimic  or  control  oligos  for  72  h,  after  which  cell  lysates 
were  harvested.  Intact  and  cleave  caspase-3  was  detected  using  an  anti-caspase-3  antibody  that  recognizes  both  the 
intact  and  the  cleaved  caspase-3.  Calnexin  levels  were  measured  as  a  loading  control.  (D-E)  Colony  formation  assay  as  a 
function  of  miR-449a  mimic.  Treatments  were  conducted  in  5  replicates.  Show  are  the  representative  images  (D)  and 
quantification  of  the  number  of  colonies  (E).  *  P  <  0.05,  comparing  to  Control. 


The  differentiation-inducing  function  of  miR-449a  is 
mediated  by  down-regulating  multiple  direct  targets  of 
miR-449a 

miRNAs  regulate  gene  expression  largely  by  decreasing  their 
target  mRNA  levels38,39  via  binding  to  specific  target  sites  in  the 


3’UTRs.  In  order  to  com¬ 
prehensively  identify  the 
targets  that  mediate  the  dif¬ 
ferentiation-inducing  func¬ 
tion  of  miR-449a,  we  first 
exploited  microarray  analy¬ 
sis  to  identify  the  predicted 
miR-449a  targets  that  are 
down-regulated  by  miR- 
449a  overexpression  in  BE 
(2)-C  cells.  Figure  3A 
shows  the  extent  of  miR- 
449a  overexpression  fol¬ 
lowing  transfection  with 
miR-449a  mimic  for  24  h. 
As  shown  in  Figure  3B, 
much  more  of  the  predicted 
miR-449a  targets  decrease 
in  expression  than  increase. 
In  contrast,  as  shown  in 
Figure  3C,  the  changes  in 
mRNAs  not  predicted  as 
miR-449a  targets  follow  a 
normal  distribution  pat¬ 
tern — roughly  equal  num¬ 
bers  of  mRNAs  not- 
predicted  as  miR-449a  tar¬ 
gets  decrease  and  increase 
in  expression.  The  empiri¬ 
cal  density  curves  in  Fig¬ 
ure  3D  further  show  the 
difference  in  the  2  distribu¬ 
tions,  and  the  difference  is 
statistically  significant 
(p  =  2.21e-30).  These 

results  indicate  that  a  sig¬ 
nificant  portion  of  the 
miR-449a  predicted  targets 
are  true  miR-449a  targets. 
The  predicted  miR-449a 
targets  that  are  downregu- 
lated  by  miR-449a  overex¬ 
pression  potentially  play 
roles  in  mediating  the  dif¬ 
ferentiation-inducing  func¬ 
tion  of  miR-449a.  In  this 
study,  we  focused  on  inves¬ 
tigating  the  top  25  pre¬ 
dicted  targets  that  are 
down-regulated  by  >40% 
by  miR-449a  mimic 
(Table  S2),  since  this  level  of  down-regulation  of  a  functional 
gene  typically  leads  to  a  dramatic,  biologically  meaningful  pheno¬ 
type.20,40"42  We  confirmed  the  downregulation  of  the  25  target 
genes  by  miR-449a  mimic  using  qPCR  (Fig.  S3).  By  exploiting  a 
set  of  siRNAs  against  the  25  genes,  we  examined  whether 
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Figure  3.  For  figure  legend,  see  page  543. 
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knocking  down  each  of  the  25  target  gene  expression  induced  cell 
differentiation  in  BE(2)-C  cells.  The  depletions  of  the  gene 
expression  by  the  corresponding  siRNAs  were  confirmed  by 
qPCR  (Fig.  S4).  As  shown  in  Figure  3E,  siRNAs  against  5  of  the 
25  target  genes  significantly  induced  neurite  outgrowth.  To 
exclude  the  possibility  that  the  neurite  outgrowth  induced  by  the 
siRNAs  is  due  to  the  off-targets  effect  of  the  synthetic  oligos,  we 
used  a  second  set  of  siRNAs  against  these  5  genes  to  validate  their 
function.  As  shown  in  Figure  3F  and  Supplementary  Figure  5, 
knockdown  of  only  3  genes,  microfibril-associated  protein  4 
(MFAP4),  Plakophilin  4  (PKP4)  and  TRNA  Splicing  Endonu¬ 
clease  15  Homolog  (TSEN15),  was  confirmed  to  induce  neurite 
outgrowth.  In  addition,  we  further  show  that  knockdown  of  the 
3  genes  dramatically  induces  expression  of  the  cell  differentiation 
markers  at  protein  levels  (Fig.  3G)  and  significantly  reduces  cell 
viability  (Fig.  3H),  further  demonstrating  that  true  cell  differen¬ 
tiation  is  induced  by  depletion  of  these  3  genes. 

To  further  characterize  the  role  of  the  3  genes  in  mediating  the 
differentiation-inducing  function  of  miR-449a,  we  investigated 
whether  overexpression  of  these  genes  block  the  differentiation 
effect  of  miR-449a  in  BE(2)-C  cells.  We  used  expression  con¬ 
structs  expressing  only  the  protein-coding  regions  of  the  genes 
(i.e.,  lacking  the  3’UTR  regions)  to  avoid  the  down-regulation  of 
the  genes  by  miR-449a.  The  protein  overexpression  levels  of  the 
3  genes  by  their  corresponding  expression  vectors  were  confirmed 
by  Western  blots,  as  shown  in  Supplementary  Figure  6.  As 
shown  Figure  31,  overexpression  of  each  of  the  genes  individually 
does  not  reduce  the  differentiation-inducing  effect  of  miR-449a, 
as  measured  by  expression  of  differentiation  markers.  However, 
co-overexpression  of  the  3  genes  dramatically  inhibits  the  effect 
of  miR-449a  on  the  expression  of  the  differentiation  markers. 
The  results  suggest  that  the  differentiation-inducing  effect  of 
miR-449a  is  mediated  by  coordinative  downregulation  of  the  3 
genes. 

Overexpression  of  miR-449a  induces  G0/G1  arrest  by 
directly  targeting  CDK6  and  LEF1 

Previous  studies  have  demonstrated  that  miR-449a  inhibits 
cell  growth  and  reduces  cell  viability  through  arresting  cells  in 
G0/G1  phase  in  other  cancer  type.28  We  therefore  examined 
whether  miR-449a  has  a  similar  function  in  neuroblastoma 


cells.  As  shown  in  Figure  4A-C,  comparing  to  the  control 
cells,  miR-449a  overexpression  leads  to  marked  increases  of 
cell  numbers  in  G0/G1  phase  and  dramatic  decreases  of  cells 
numbers  in  S  and  G2/M  phases,  suggesting  that  miR-449a 
blocks  cell  cycle  progression  through  arresting  cells  in  G0/G1 
phase.  We  further  show  that  the  effect  of  miR-449a  on  cell 
cycle  distribution  is  comparable  to  that  induced  by  RA.  It  has 
been  known  that  RA  treatment  of  neuroblastoma  cells  lead  to 
terminal  cell  differentiation,  which  is  accompanied  with  cell 
cycle  exit.43  We  therefore  examined  whether  miR-449a  overex¬ 
pression  has  similar  effect  in  neuroblastoma  cells  by  examining 
the  BrDU  incorporation  into  DNA  in  the  nucleus  of  individ¬ 
ual  cells  using  immunofluorescence  staining.  As  shown  in  Fig¬ 
ure  4D-G,  comparing  to  the  control  cells,  more  cells  treated 
with  miR-449a  mimic  or  RA  are  negative  for  BrDU  staining, 
which  indicates  that  these  cells  failed  to  enter  S-phase.  These 
results  indicate  that  induction  of  cell  cycle  exit  is  achieved  by 
miR-449a  overexpression  in  neuroblastoma  cells,  and  the 
effect  of  miR-449a  in  inducing  cell  cycle  exit  is  comparable  to 
RA. 

We  further  intended  to  identify  the  targets  that  mediate 
the  effect  of  miR-449a  on  cell  cycle  arrest.  Among  the  25  tar¬ 
gets  that  were  downregulated  greater  than  40%  by  miR-449a 
mimic,  7  genes,  including  CCNE2,  CDC25A,  CDK6,  LEF1, 
STMN1,  CASP2  and  HDAC1,  are  identified  by  Ingenuity 
Pathway  Analysis  (IPA,  www.ingenuity.com)  analysis  as 
involved  in  regulating  cell  cycle  progression.  We  first  exam¬ 
ined  whether  individual  knockdown  of  these  7  genes  can 
exert  an  effect  on  cell  viability.  Surprisingly,  as  shown  in  Fig¬ 
ure  5A-B,  only  the  knockdown  of  CDK6  and  LEF1  signifi¬ 
cantly  reduces  cell  viability  in  BE(2)-C  and  SKNBE  cells. 
The  depletions  of  the  gene  expression  by  the  corresponding 
siRNAs  in  BE(2)-C  cells  were  confirmed  by  qPCR  (Fig.  S4). 
This  effect  of  CDK6  and  LEF1  knockdown  on  cell  survival 
was  further  confirmed  using  a  2nd  set  of  siRNAs  in  BE(2)-C 
and  SKNBE  cells  (Fig.  5C-D),  and  was  validated  in  KELLY 
cells  using  2  sets  of  siRNAs  (Fig.  5E).  Further  examination 
of  cell  cycle  distributions  indicated  that  knocking  down 
CDK6  and  LEF1  indeed  leads  to  increased  accumulation  of 
cells  in  G0/G1  phase  as  compared  with  control  oligo 
(Fig.  5F-H)  in  all  the  tree  cell  lines  BE(2)-C,  SKNBE  and 


Figure  3  (See  previous  page).  The  systematic  identification  of  miR-449a  targets  that  regulate  neuroblastoma  cell  differentiation.  (A)  Expression 
level  of  miR-449a  in  BE(2)-C  cells  after  transfected  with  miR-449a  (25  nM)  mimic  for  24  h,  as  measured  by  qPCR.  (B-C)  Histogram  of  changes  in  mRNA 
expression  levels  of  genes  (B)  predicted  as  miR-449a  targets  and  (C)  not  predicated  as  miR-449a  targets  induced  by  miR-449a  overexpression.  The  mRNA 
expression  levels  were  determined  by  gene  expression  microarray.  Shown  is  the  relative  frequency  of  observed  log2  fold  changes  for  the  predicted  miR- 
449a  targets.  (D)  The  empirical  density  curves  generated  based  on  the  above  frequency  distributions  by  Gaussian  Kernel  Smoothing.  (E)  Effect  of  knock¬ 
ing  down  the  top  25  target  genes  on  neurite  outgrowth  in  BE(2)-C  cells.  Cells  were  transfected  with  25  nM  of  the  indicated  siRNAs  or  negative  control  oli¬ 
gos  (Control)  in  triplicated  for  4  d.  Neurite  outgrowth  was  measured  as  above.  *,  P  <  0.05,  comparing  to  Control.  (F)  Validation  of  the  effect  of  knocking 
down  the  5  candidate  targets  on  neurite  outgrowth  using  the  2nd  set  siRNAs.  BE(2)-C  cells  were  transfected  and  neurite  outgrowth  were  measured  as 
above.  *  P  <  0.05,  comparing  to  control.  (G)  Effect  of  gene  knockdown  on  expression  of  differentiation  markers  in  BE(2)-C  cells.  Cells  were  transfected 
with  the  indicated  siRNAs  or  control  oligo  (25  nM).  Protein  levels  were  measured  by  Western  blot  as  above.  (H)  Effect  of  gene  knockdown  on  cell  viability 
in  BE(2)-C  cells.  Cells  were  transfected  with  the  indicated  siRNAs  or  control  oligo  (25  nM)  in  triplicates  in  96-well  plates.  After  4  d,  cell  viability  was  mea¬ 
sured  as  above.  Shown  are  the  relative  cell  viabilities,  presented  as  averages  of  the  3  replicates,  normalized  to  the  control  oligos.  *,  P  <  0.05,  comparing 
to  the  control.  (I)  Effect  of  gene  overexpression  on  expression  of  differentiation  markers.  BE(2)-C  cells  were  co-transfected  with  the  indicated  expression 
vectors  and  oligos.  After  4  d,  protein  levels  of  expression  markers  were  measured  as  above.  *,  P  <  0.05,  comparing  to  Control. 
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Figure  4.  Effect  of  miR-449a  on  cell  cycle  progression  in  neuroblastoma  cells.  (A-C)  Effect  of  miR-449a  overexpression  on  cell  cycle  distribution  in  (A) 
BE(2)-C,  (B)  SKNBE  and  (C)  KELLY  cells.  Cells  were  transfected  with  25  nM  miR-449a  mimic  or  control  oligo,  or  treated  with  5  |jlM  RA.  Mock  transfected 
cells  were  used  as  an  additional  negative  control.  After  3  d,  cells  were  collected  and  stained  with  propidium  iodide  for  cell  cycle  analysis.  The  fraction  of 
cells  in  G0/G1,  S  and  G2  phases  was  quantified  using  the  Jean-Jett-Fox  model.  Similar  results  were  obtained  from  3  independent  experiments.  (D-G) 
Effect  of  miR-449a  and  RA  on  DNA  synthesis  as  measured  by  BrDU  incorporation  in  BE(2)-C  (D),  SKNBE  (E)  and  KELLY  (F)  cells.  Cells  were  transfected  with 
miR-449a  mimic  or  control  oligos  (25  nM),  or  treated  with  RA  (5  |jlM)  for  4  d.  Cells  were  incubated  with  BrdU  during  the  final  2  h  of  culture.  Cells  were 
then  stained  for  BrDU  incorporation  (red)  as  described  in  the  Materials  and  Methods.  Cell  nuclei  were  identified  by  nuclear  staining  with  DAPI  (blue).  (G) 
Quantification  of  percentage  of  cells  positive  for  BrDU  staining.  *,  P  <  0.05,  comparing  to  control. 
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Figure  5.  Function  of  the  predicted  miR-449a  cell  cycle-regulating  targets  in  regulating  cell  cycle  progression  and  cell  survival.  (A-E),  Effect  of 
knocking  down  the  predicted  miR-449a  targets  on  cell  viability  in  BE(2)-C,  SKNBE  and  KELLY  cells  using  2  sets  of  siRNAs  for  the  indicated  genes.  Cells 
were  transfected  with  25  nM  of  the  indicated  siRNAs  or  control  oligo  (Control).  After  4  d,  cell  viability  was  measured  as  above.  *,  P  <  0.05,  comparing  to 
Control.  (F-H),  Effect  of  knocking  down  CDK6  and  LEF1  on  cell  cycle  distribution  in  (F)  BE(2)-C,  (G)  SKNBE  and  (H)  KELLY  cells.  Cells  were  transfected  with 
25  nM  of  the  indicated  siRNAs  or  control  oligo  (Control).  Mock-transfected  cells  were  used  as  an  additional  negative  control.  After  3  d,  cell  cycle  distribu¬ 
tions  were  measured  as  above.  (I-K)  Effect  of  knocking  down  CDK6  and  LEF1  on  BrDU  incorporation  in  (I)  BE(2)-C,  (J)  SKNBE  and  (K)  KELLY  cells.  Cells 
were  transfected  as  above,  and  the  level  of  BrDU  incorporation  was  measured  as  above.  (L)  Western  blot  analysis  of  Caspase-3  activation.  BE(2)-C  cells 
were  transfected  with  25  nM  of  the  indicated  oligos  for  3  d,  and  caspase-3  levels  were  measured  as  above.  (M)  Effect  of  gene  overexpression  on  cell  cycle 
distribution.  BE(2)-C  cells  were  transfected  with  the  indicated  expression  vectors  and  oligos.  After  3  d,  cell  cycle  distributions  were  analyzed  as  above. 
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Figure  6.  Validation  of  the  predicted  direct  targets  of  miR-449a.  (A)  The  predicted  interactions  between  miR- 
449a  and  the  target  sites  in  the  3'UTRs  of  the  indicated  mRNAs.  The  seed  sequences  are  underlined.  (B)  Validation 
of  the  target  sites  of  miR-449a  by  luciferase  reporter  assay.  BE(2)-C  cells  were  co-transfected  with  the  indicated  vec¬ 
tors  and  miR-449a  mimic.  After  3  d  of  transfection,  cells  were  lysed  and  luciferase  activity  was  measured.  Shown  are 
normalized  luciferase  activities  with  the  luciferase  activity  associated  with  wildtype  (WT)  3'UTRs  normalized  to  those 
associated  with  the  corresponding  mutant  3'UTRs.  *,  P  <  0.05,  comparing  to  mutant  3'UTRs.  (C)  Validation  of  the 
direct  interactions  of  miR-449a  with  the  target  genes  by  Bio-miR-449a  pulldown  assay.  BE(2)-C  cells  were  transfected 
with  the  25  nM  of  the  Bio-miR-449a  or  control  oligo  (Bio-Control).  The  pulldown  assay  was  performed  following  the 
protocol  described  in  the  Materials  and  Methods.  Shown  are  fold  of  enrichment  of  the  indicated  mRNAs  by  Bio-miR- 
449a,  calculated  by  normalizing  the  levels  of  mRNAs  associated  with  Bio-miR-449a  to  those  associated  with  the  Bio¬ 
control.  *,  P  <  0.05,  comparing  to  the  Bio-control.  (D)  BIRC5  and  CDT1  mRNA  levels  in  the  Biotin  pulldown  products 
were  measured  by  qPCR.  *,  P  <  0.05,  comparing  to  the  Bio-control. 


CDK6  31  UTR 

5'  UACUUUCUGCCACACACTJGCCU 
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KELLY.  Correspondingly, 

CDK6  and  LEL1  knock¬ 
down  leads  to  reduced  cell 
proliferation,  as  measured 
by  BrDU  incorporation 
assay  (Fig.  5I-K).  In  addi¬ 
tion,  we  show  that  CDK6 
and  LEF 1  knockdown 
also  leads  to  the  activation 
of  the  apoptotic  pathway 
as  indicated  by  the 
increased  levels  of  cleaved 
capase  3  (Fig.  5L),  consis¬ 
tent  with  the  effect  of 
miR-449a  overexpression. 

We  further  examined  the 
role  of  CDK6  and  LEF1 
in  mediating  the  cell 
cycle-regulating  function  of 
miR-449a  by  overexpressing 
these  2  genes  in  BE(2)-C 
cells.  The  protein  overex¬ 
pression  levels  of  the  2  genes 
by  their  corresponding 
expression  vectors  were  con¬ 
firmed  by  Western  blots 
(Fig.  S6).  As  shown  in  Fig¬ 
ure  5M,  overexpression  of 
CDK6  and  LEF1  individu¬ 
ally  decrease  the  effect  of 
miR-449a  on  G0/G1  arrest, 
and  co-overexpression  of 
CDK6  and  LEF1  shows 
enhanced  effect  comparing 
to  individual  gene  overex¬ 
pression.  These  results  indi¬ 
cates  the  function  of  miR- 
449a  in  regulating  cell  cycle 
distribution  is  mediated  by 
coordinately  down-regulat¬ 
ing  expression  of  CDK6 
and  LEF  1 . 

Overall,  our  results  sug¬ 
gest  that,  among  the  7  tar¬ 
get  genes  involved  in  cell 
cycle  progression,  CDK6 
and  LEF1  play  dominant 
roles  in  mediating  the  effect 
of  miR-449a  on  cell  cycle 
progression. 

miR-449a  directly  targets  the  3'  UTRs  of  MFAP4,  PKP4, 
TSEN15,  CDK6  and  LEF1 

The  above  results  identified  3  direct  targets  of  miR-449a  that 
mediate  its  effect  on  cell  differentiation  and  2  targets  that  mediate 
its  effect  on  cell  cycle  arrest.  The  predicted  target  sites  of  miR- 


449a  in  the  mRNA  3/UTRs  are  shown  in  Figure  6A.  We  used  a 
luciferase  reporter  assay  to  validate  the  predicted  target  sites.  As 
shown  in  Figure  6B,  miR-449a  over-expression  significantly 
decreased  the  luciferase  activity  in  cells  expressing  the  wildtype 
3rUTRs  as  compared  with  mutated  3rUTRs  for  all  of  the  5  target 
genes.  For  each  of  the  targets,  the  luciferase  activity  for  wildtype 
3rUTRs  were  decreased  by  >25  %  comparing  to  the  mutated 
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Figure  7.  Correlation  of  tumor  mRNA  expression  levels  with  survival  of  neuroblastoma  patients  in  the  Vertseeg  patient  cohort.  For  each  of  the 
indicated  gene,  the  88  neuroblastoma  patients  were  divided  into  high  and  low  groups  as  described  in  the  Materials  and  Methods.  Shown  are  the  mRNA 
expression  levels  (A,  C,  E,  G,  I)  and  the  Kaplan-Meier  overall  survival  curves  (B,  D,  F,  H,  J)  in  the  low  and  high  expression  groups  for  each  of  the  indicated 
genes.  The  mRNA  expression  levels  were  presented  as  average  levels  of  the  normalized  signal  intensities  in  each  group. 


3/UTRs.  This  magnitude  of  repression  is  consistent  with  that 
observed  in  previous  studies  of  miRNA  regulation  using  lucifer- 
ase  reporter  vectors.44,45  These  results  demonstrate  that  the  pre¬ 
dicted  target  sites  are  true  target  sites  of  miR-449a.  In  order  to 
further  demonstrate  the  physical  interactions  of  these  mRNAs 
with  miR-449a,  we  performed  biotinylated-miRNA  pulldown 
assay46  in  BE(2)-C  cells  transfected  with  biotinylated  miR-449a 
mimic  (Bio-miR-449a).  As  shown  in  Figure  6C,  for  all  the  5  tar¬ 
get  genes  examined,  their  mRNA  levels  are  significantly  enriched 


by  Bi-miR-449a  pull-down  relative  to  the  biotinylated  control 
oligo  (Bio-control).  In  contrast,  the  mRNA  levels  of  BIRC5 
(Baculoviral  IAP  Repeat  Containing  5)  and  CDTl  (Chromatin 
Licensing  and  DNA  Replication  Factor  1),  which  are  not  targets 
of  miR-449a,  are  not  enriched  by  Bio-miR-449a  pulldown 
(Fig.  6D),  further  supporting  that  Bio-miR-449a  pulldown 
product  is  specifically  enriched  with  mRNAs  that  are  miR-449a 
targets.  Overall,  our  studies  demonstrate  the  physical  interactions 
of  the  5  target  mRNAs  with  miR-449a. 
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High  tumor  expression  levels  of  MFAP4,  PKP4  and 
TSEN15,  the  differentiation-regulating  targets  of  miR-449a, 
are  correlated  with  poor  survival  of  neuroblastoma  patients 

The  above  results  suggest  that  miR-449a  functions  as  a  tumor 
suppressor  in  neuroblastoma  through  at  least  2  mechanisms — 
inducing  cell  differentiation  and  arresting  cells  in  G0/G1  phase. 
In  order  to  further  understand  the  clinical  relevance  of  the  miR- 
449a-mediated  tumor  suppressive  pathway  in  neuroblastoma 
development,  we  assessed  the  correlation  of  the  expression  of  the 
identified  target  genes  in  neuroblastoma  tumor  specimens  with 
patient  survival  in  2  published  neuroblastoma  patient  cohorts: 
the  Versteeg  and  the  Kocak  cohorts.47,48  As  shown  in  Fig¬ 
ure  7A-F  and  Figure  8A-F,  Kaplan-Meier  survival  analysis  indi¬ 
cates  that,  in  both  patient  cohorts,  patients  with  high  tumor 
MFAP4,  PKP4  or  TSEN15  mRNA  levels  have  significantly 
lower  overall  survival  rates  comparing  to  patients  with  low 
expression  levels  of  these  genes.  However,  as  shown  in  Fig¬ 
ure  7G-J  in  the  Versteeg  cohort,  tumor  expression  levels  of 
CDK6  and  LEF1  are  not  significantly  correlated  with  patient  sur¬ 
vival.  We  observe  similar  results  for  CDK6  in  the  Kocak  cohort 
(Fig.  8G-H).  Paradoxically,  in  the  Kocak  cohort,  the  low  levels 
of  LEF1  in  tumors  are  correlated  with  poor  survival  (Fig.  8I-J), 
which  contradicts  the  tumor  suppressive  function  of  miR-449a. 
Overall,  our  comprehensive  analyses  suggest  that  expression  levels 
of  the  identified  differentiation-regulating  targets  of  miR-449a, 
but  not  the  cell  cycle-regulating  targets,  are  clinically  relevant  in 
determining  the  prognosis  of  the  neuroblastoma  patients. 

Discussion 

Our  study  identifies  a  novel  tumor  suppressive  mechanism  in 
neuroblastoma  mediated  by  miR-449a.  Our  results  clearly  indi¬ 
cate  that  the  function  of  miR-449a  in  reducing  neuroblastoma 
cell  survival  involves  3  molcular  mechanisms:  inducing  cell  cycle 
arrest  and  cell  differentiation,  inducing  cell  apoptosis,  and  reduc¬ 
ing  cell  proliferation.  Previous  studies  have  demonstrated  that 
miR-449a  functions  as  a  tumor  suppressor  through  blocking  cell 
cycle  progression  in  several  types  of  cancers,  including  lung  can¬ 
cer,  prostate  cancer,  gastric  cancer  and  retinoblastoma.26"29,32 
However,  whether  miR-449a  regulates  cell  cycle  progression  in 
neuroblastoma  has  not  been  investigated.  Consistent  with  previ¬ 
ous  studies  in  other  cancer  types,  we  show  that  miR-449a  also 
induces  G0/G1  phase  arrest  in  neuroblastoma  cells.  Aside  from 
its  function  in  regulating  cell  cycle,  we  show  here  that  miR-449a 
functions  as  a  tumor  suppressor  through  an  additional  mecha¬ 
nism — inducing  neuroblastoma  cell  differentiation.  In  addition, 
we  identified  novel  target  genes  of  miR-449a  that  mediate  its 
function  in  regulating  cell  differentiation  and  cell  cycle  progres¬ 
sion.  It  is  possible  that  the  reduced  cell  proliferation  and  induced 
cell  apoptosis  by  miR-449a  overexpression  are  just  the  subse¬ 
quent  consequences  of  cell  differentiation  and  cell  cycle  arrest. 
However,  it  is  possible  that  miR-449a  also  targets  genes  in  the 
cell  proliferation  and  apoptosis  pathways  that  are  independent  of 
the  cell  cycle/differentiation-regulating  pathways.  This  is  cer¬ 
tainly  interesting  and  warrants  further  investigations  in  the 


Figure  8.  Correlation  of  tumor  mRNA  expression  levels  with  survival 
of  neuroblastoma  patients  in  the  Kocak  patient  cohort.  The  476 

patients  were  divided  into  high  and  low  groups  as  above.  Shown  are  the 
mRNA  expression  levels  (A,  C,  E,  G,  I)  and  the  Kaplan-Meier  overall  sur¬ 
vival  curves  (B,  D,  F,  H,  J)  in  the  low  and  high  expression  groups  for 
each  of  the  indicated  genes.  The  mRNA  expression  levels  were  presented 
as  average  levels  of  the  normalized  signal  intensities  in  each  group. 


future.  Overall,  our  study  characterizes  for  the  first  time  the 
underlying  the  molecular  mechanisms  that  mediate  the  tumor 
suppressive  function  of  miR-449a  in  neuroblastoma  cells,  as  illus¬ 
trated  in  Figure  9.  Our  in  vitro  investigations  certainly  suggest 
that  this  tumor  suppressive  mechanism  mediated  miR-449a  may 
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Figure  9.  The  mechanisms  of  the  miR-449a  tumor  suppressive  function  in  neuroblastoma 
cells.  Our  results  support  that  miR-449a  functions  as  a  tumor  suppressor  in  neuroblastoma 
through  2  molcular  pathways:  The  cell  differentiation  pathway,  which  is  mediated  by  down-regu¬ 
lating  its  direct  targets  MFAP4,  PKP4  and  TSEN15;  The  cell  cycle  pathway,  which  is  mediated  by 
down-regulating  its  direct  targets  CDK6  and  LEFI.  The  induced  cell  differentiation  and  cell  cycle 
arrest  will  result  in  reduced  cell  proliferation  and  increased  cell  apoptosis,  which  will  lead  to 
tumor  growth  arrest.  However,  it  is  possible  that  miR-449a  also  targets  genes  in  the  cell  prolifera¬ 
tion  and  apoptosis  pathways  that  are  independent  of  the  cell  cycle/differentiation-regulating 
pathways — this  potential  mechanism  needs  to  be  further  investigated  in  the  future. 


play  a  critical  role  in  neuroblastoma  tumorigenicity.  However, 
future  in  vivo  investigations  in  neuroblastoma  animal  models  are 
clearly  needed  in  order  to  fully  define  the  role  of  this  miR-449a- 
mediated  molecular  pathway  in  determining  neuroblastoma  initi¬ 
ation  and  progression. 

In  our  investigation  of  multiple  neuroblastoma  cell  lines  with 
different  genetic  backgrounds,  we  found  that  although  miR-449a 
shows  a  general  effect  on  cell  differentiation,  the  sensitivity  of  dif¬ 
ferent  cell  lines  to  miR-449a-induced  differentiation  varies 
greatly — it  only  induce  cell  differentiation  to  a  minor  extent  in 
some  cell  lines  (LAN6  and  KELLY  cells).  Our  previous  work 
shows  that  another  miRNA,  miR-506-3p  shows  a  strong  differ¬ 
entiation-inducing  effect  in  KELLY  cells.20  The  sensitive 
response  of  KELLY  cells  to  the  miR-506-3p-mediated  differenti¬ 
ation-inducing  pathway  indicates  that  this  cell  line  does  not  lack 
intrinsic  differentiation  capacity,  suggesting  that  its  lack  of 
response  to  miR-449a  is  a  miR-449a-specific  effect.  The  observed 
distinct  response  of  a  particular  neuroblastoma  cell  line  to  differ¬ 
ent  differentiation-inducing  miRNAs  is  certainly  interesting  and 
potentially  has  great  clinical  significance.  It  not  only  indicates 
that  different  differentiation-inducing  miRNAs  functions 
through  distinct  molecular  pathways,  but  also  suggest  that  neuro¬ 
blastoma  cell  differentiation  can  be  induced  through  different 
molecular  pathways  by  differentiation  agents  that  act  through 
distinct  molecular  mechanisms.  Lor  neuroblastoma  that  is  not 
responsive  to  one  type  of  differentiation  agents,  it  is  possible  to 
use  differentiation  agents  of  different  mechanisms  to  successfully 
treat  it. 

Interestingly,  despite  their  distinct  responses  to  the  differentia¬ 
tion  pathway,  all  the  cell  lines  respond  well  to  miR-449a-induced 
cell  growth  arrest.  These  results  suggest  that,  when  one  tumor 
suppressive  pathway  mediated  by  miR-449a  is  inactive,  the  other 
can  overcome  the  silencing  of  this  pathway  to  achieve  its  strong 
tumor  suppressive  function — this  dual  tumor  suppressive  mecha¬ 
nisms  make  miR-449a-based  therapeutics  a  promising  approach 


to  treat  neuroblastoma  with  a  wide  spec¬ 
trum  of  genetic  backgrounds. 

Our  study  identifies  novel  targets  of 
miR-449a  that  mediate  its  tumor  suppres¬ 
sive  function.  It  is  well  known  that  each 
miRNA  has  hundreds  of  direct  targets,49 
and  each  biological  phenotype  induced  by  a 
specific  miRNA  is  likely  achieved  by  con- 
cordantly  down-regulating  multiple  tar¬ 
gets.50”52  Multiple  miR-449a  targets  that 
are  involved  in  mediating  its  tumor  suppres¬ 
sor  function  have  been  identified  and  vali¬ 
dated  in  other  cancer  types.  However,  our 
investigation  shows  that  only  8  of  the  20 
targets  validated  in  other  cancer  cell  types 
are  down-regulated  in  neuroblastoma  cells 
(Table  S3).  This  suggests  that  the  function 
of  miR-449a  in  regulating  its  target  expres¬ 
sion  is  cell-context  specific.  More  impor¬ 
tantly,  we  identified  novel  targets  of  miR- 
449a,  including  MLAP4,  PKP4  and 
TSEN15,  that  play  roles  in  regulating  neuroblastoma  cell  differ¬ 
entiation.  The  molecular  functions  of  these  3  genes  have  been 
defined  in  previous  studies.53”57  However,  the  role  of  these  genes 
in  tumorigenesis  as  well  as  in  regulating  cell  differentiation  has 
not  been  intensively  explored  previously — we  are  the  first  to 
identify  their  potential  oncogenic  functions  in  neuroblastoma  by 
participating  the  cell  differentiation  process.  More  remarkably, 
we  found  that  all  the  3  differentiation-regulating  targets  of  miR- 
449a,  MLAP4,  PKP4  and  TSEN15,  are  significantly  correlated 
with  neuroblastoma  patient  survival,  with  high  mRNA  expres¬ 
sion  levels  of  these  genes  in  neuroblastoma  tumor  specimens  cor¬ 
related  with  poor  patient  survival.  These  findings  strongly 
suggest  their  critical  role  in  determining  the  prognosis  of  neuro¬ 
blastoma  patients.  Lurther  studies  are  warranted  to  characterize 
the  molecular  mechanisms  by  which  these  genes  regulate  neuro¬ 
blastoma  cell  differentiation,  to  elucidate  how  they  fit  into  the 
known  signaling  network  that  control  neuroblastoma  cell  differ¬ 
entiation,  and  to  investigate  whether  specific  differentiation 
agents  targeting  these  genes  can  be  developed. 

Of  note,  in  this  study,  we  found  that  the  expression  levels  of 
CDK6  in  neuroblastoma  are  not  significantly  correlated  with 
patient  survival,  although  we  observed  its  role  in  regulating  cell 
cycle  progression  in  vitro  (Fig.  5)  and  their  oncogenic  functions 
have  been  well  demonstrated  previously.58,59  One  possible  expla¬ 
nation  is  that  CDK6  may  need  to  coordinate  with  additional 
oncogenic  pathways  in  neuroblastoma  cells  in  order  to  reach  a 
clinical  significant  impact  on  patient  survival.  This  is  certainly  an 
interesting  possibility  and  worth  to  be  further  pursued.  The 
oncogenic  function  of  LEFI  has  also  been  reported  previously.60" 
62  However,  we  observed  contradictive  results  in  2  neuroblas¬ 
toma  patient  cohorts,  which  leaves  the  association  of  LEFI 
expression  with  neuroblastoma  patient  prognosis  undefined  in 
our  study.  Future  studies  are  certainly  needed  to  clearly  define 
the  function  of  LEFI  in  neuroblastoma  development.  In  addi¬ 
tion,  the  role  of  CDK6  and  LEFI  in  regulating  cell 
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differentiation  has  also  been  indicated  in  previous  studies.63,64 
However,  we  failed  to  observe  the  function  of  these  2  genes  in 
regulating  neuroblastoma  cell  differentiation.  These  results  fur¬ 
ther  demonstrate  the  complexity  and  the  cell-context  specificity 
of  the  cell  differentiation  pathways — it  is  reasonable  to  believe 
that  there  are  alternative,  cell  type-specific  signaling  pathways  in 
neuroblastoma  cells  to  control  the  cell  differentiation  process 
that  are  normally  controlled  by  CDK6  and  LEF1  in  other  cell 
types,  which  makes  CDK6  and  LEF1  not  essential  to  determine 
the  differentiation  fate  of  neuroblastoma  cells. 

Another  finding  in  our  study  is  that  co-overexpression  of  the  3 
differentiation-regulating  targets  MFAP4,  PKP4  and  TSEN15 
only  partially  inhibited  the  differentiation-inducing  effect  of 
miR-449a.  Likewise,  co-overexpression  of  CDK6  and  LEF 1  only 
partially  inhibited  the  effect  of  miR-449a  on  cell  cycle  distribu¬ 
tion.  These  results  suggest  that  there  are  additional  targets  play¬ 
ing  an  important  role  in  its  differentiation-inducing  and  cell 
cycle-regulating  functions.  In  this  study,  we  only  investigated  the 
genes  that  are  predicted  as  miR-449a  targets  using  the  canonical 
miRNA  target  prediction  approach,  which  is  based  on  the  inter¬ 
actions  of  seed  region,  a  6-8  nucleotides  at  the  5'  end  of  the 
miRNA,  with  the  target  sites  in  the  mRNA  3'  UTR, 65-67  and  we 
only  investigated  the  predicted  targets  genes  that  are  downregu- 
lated  by  >40%  at  mRNA  levels  by  miR-449a.  Although  miR- 
NAs  down-regulate  the  majority  of  their  targets  at  mRNA  levels, 
it  is  known  miRNAs  can  also  regulate  its  target  gene  expression 
through  translational  repression,  leading  to  decreased  protein 
expression  of  the  target  genes  without  affecting  the  mRNA  lev¬ 
els.68  In  addition,  it  has  been  demonstrated  that  miRNA  may 
also  target  genes  through  non-canonical  interactions.69-72  It  is 
highly  possible  that  there  are  additional  miR-449a  targets  that 
not  investigated  in  this  study  playing  important  roles  in  mediat¬ 
ing  its  differentiation-inducing  and  cell  cycle-regulating  func¬ 
tions.  Future  studies  are  certainly  warranted  to  investigate  the 
role  of  these  targets  in  order  to  completely  elucidate  the  tumor 
suppressive  mechanisms  of  miR-449a  in  neuroblastoma  cells. 

Due  to  the  limited  resources  for  investigating  the  miR-449a 
expression  in  neuroblastoma  tumor  specimens,  we  were  unable  to 
examine  the  miR-449a  endogenous  expression  levels  in  neuroblas¬ 
toma  primary  tumors  and  therefore  unable  to  characterize  the  cor¬ 
relation  of  miR-449a  levels  with  neuroblastoma  patient  survival  in 
the  present  study.  However,  we  show  that  the  endogenous  expres¬ 
sion  of  miR-449a  in  differentiated  neuroblastoma  cells  is  signifi¬ 
cantly  increased  comparing  to  the  undifferentiated  cells.  These 
results  strongly  suggest  miR-449a  is  endogenously  expressed  in 
normal  differentiated  neural  crest  cells,  but  its  expression  is  signifi¬ 
cantly  suppressed  in  undifferentiated  neuroblastoma  cells — this 
raises  the  possibility  to  activate  its  endogenous  expression  as  a  pos¬ 
sible  approach  to  differentiation  therapy.  The  mechanisms  that 
control  its  endogenous  expression  certainly  warrant  further  investi¬ 
gation  in  order  to  develop  therapeutic  strategies  to  activate  its 
endogenous  differentiation-inducing  function. 

To  summarize,  we  have  identified  a  novel  miR-449a-mediated 
tumor  suppressive  mechanism  in  neuroblastoma.  The  clinically 
observed  significant  correlation  of  the  differentiation-regulating 
targets  of  miR-449a  with  neuroblastoma  patient  survival  suggests 


the  potential  clinical  significance  of  the  dysregulation  of  miR- 
449a  expression  in  neuroblastoma  tumorigenesis  —  loss  of  miR- 
449a  expression  in  neuroblastoma  cells  may  play  a  critical  role  in 
neuroblastoma  development  by  inducing  aberrant  expression  of 
MFAP4,  PKP4  and  TSEN15.  Further  studies  are  definitely  war¬ 
ranted  to  investigate  the  correlation  of  tumor  miR-449a  expres¬ 
sion  with  patient  prognosis. 

Materials  and  Methods 

Materials 

miR-449a  mimic  and  negative  control  oligos  were  obtained 
from  Thermo  Fisher  Scientific.  miR-449a  precursors  were  pur¬ 
chased  from  Ambion.  Rabbit  anti-GAP43,  anti-NSE  and 
anti-|3III-tubulin  antibodies  were  obtained  from  Abeam.  Rabbit 
anti-calnexin,  and  goat  anti-rabbit  secondary  antibody  conju¬ 
gated  with  horseradish  peroxidase  (HRP)  were  purchased  from 
Santa  Cruz  Biotechnology,  Inc.  Rabbit  anti-human  caspase-3, 
CDK6  and  LEF1  antibodies  were  from  Cell  Signaling  Technol¬ 
ogy.  Rabbit  anti-human  MFAP4  and  PKP4  antibodies  were 
from  Thermo  Fisher  Scientific,  Inc.  Rabbit  anti-human  TSEN15 
antibody  was  from  Aviva  Systems  Biology.  Bio-miR-449a  and 
Bio-control  oligo  were  purchased  from  Sigma.  siRNAs  and  the 
negative  control  oligos  were  purchased  from  Dharmacon  and 
Sigma,  respectively.  pCMV-CDK6  expression  vector  was  a  gift 
from  Sander  van  den  Heuvel  (Addgene  plasmid  #1868). 73 

Cell  lines 

BE(2)-C,  SKNBE  and  BE(2)-M17  cells  were  purchased  from 
the  American  Type  Culture  Collection  (ATCC).  LAN6  cells 
were  obtained  from  Children’s  Oncology  Group.  KELLY  cells 
were  obtained  from  the  cell  line  repository  at  the  Greehey  Child¬ 
ren’s  Cancer  Research  Institute,  University  of  Texas  Health  Sci¬ 
ence  Center  at  San  Antonio.  Cells  were  grown  in  DMEM/F12 
supplemented  with  10%  fetal  bovine  serum. 

Detection  and  quantification  of  neurite  outgrowth 

Cells  were  plated  and  treated  in  96-well  plates.  For  measuring 
neurite  outgrowth,  cells  were  placed  into  IncuCyte  ZOOM  Live 
Cell  Imaging  System  (Essen  BioScience),  and  images  were  taken 
under  20X  magnification.  For  measuring  neurite  outgrowth  in  a 
time-dependent  manner,  cell  images  were  taken  every  6  hours. 
The  neurite  lengths  associated  with  each  treatment  were  calcu¬ 
lated  using  the  neurite  definition  determined  for  each  specific  cell 
line  using  the  NeuroTrack  system  (Essen  BioScience). 

Cell  cycle  analysis 

Cells  were  detached  by  scrapping,  and  cell  suspension  was 
pipetted  for  10  times  to  separate  the  cell  clumps  into  single  cells. 
Cells  were  then  collected  by  centrifuging  at  1,000  rpm  for 
5  min.  The  cells  were  washed  once  with  IX  PBS.  To  fix  cells, 
ice-cold  IX  PBS  containing  1  mM  EDTA  and  85%  ethanol  was 
added  to  cells  drop-wise  while  vortexing,  in  order  to  prevent  cell 
aggregation.  After  1  hour,  the  cells  were  harvested  by  centrifuga¬ 
tion  at  1,400  rpm  for  5  minutes  at  4°C,  and  cells  were  re-sus- 
pended  in  IX  PBS,  and  treated  with  50  p,g/ml  Propidium  iodide 
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and  100  (xg/ml  RNase  A  for  30  minutes  at  37° C.  Cells  were  then 
passed  through  271  /2G  needle  for  5  times  and  filtered  through 
35  pM  strainers.  Cell  cycle  data  were  collected  on  a  Cytomics 
FC  500  flow  cytometer  (Beckman  Coulter,  Brea,  CA).  Data  were 
analyzed  using  Flowjo  version  7.6.5  (TreeStar,  Ashland,  OR). 

Western  blots 

Cell  lysates  were  prepared  using  RIPA  buffer.  Protein  concen¬ 
tration  was  determined  using  the  Pierce  BCA  assay  (Thermo 
Fisher  Scientific).  For  electrophoresis,  equal  amounts  of  cell 
lysates  were  resolved  by  SDS-PAGE  and  transferred  to  Immun- 
Blot  PVDF  membranes  (Bio-Rad  Laboratories).  Membranes 
were  blocked  and  probed  with  primary  antibodies  to  specific  pro¬ 
teins.  Bound  antibodies  were  detected  with  secondary  antibodies 
conjugated  with  horseradish  peroxidase  (HRP)  and  visualized 
using  SuperSignal  West  Pico  Chemiluminescent  Substrate 
(Thermo  Fisher  Scientific). 

Cell  viability  assay 

Cell  viability  was  measured  as  previously  described.41  Briefly, 
cells  were  plated  in  96-well  format  and  treated  as  specified.  After 
culturing  for  4-5  d,  cell  viability  was  determined  using  the  Cell- 
Titer-Glo  Luminescent  Cell  Viability  Assay  (Promega). 

BrdU  incorporation  assay 

Cells  were  seeded  in  96-well  plates  and  treated  as  indicated. 
After  culturing  for  3  d,  cells  were  incubated  with  10  |jlM  BrdU 
(Roche,  Indianapolis,  IN,  USA),  and  then  fixed  with  LixDenat 
solution.  Anti-BrdU-POD  was  added  to  fixed  cells,  and  substrate 
solution  containing  tetramethylbenzidine  was  used  to  detect  and 
quantify  the  amount  of  BrdU  incorporation. 

Immunofluorescence  detection  of  DNA  synthesis  by  BrdU 
labeling 

We  used  a  protocol  described  by  Feliciano  DM.43  Briefly  1  x 
104  BE(2)-C  cells  were  plated  on  coverslips  in  12-well  plates,  and 
were  transfected  with  miR-449a  mimic  (25  nM)  or  control  oli- 
gos  (25  nM),  or  treated  with  RA  (5  fxM)  for  4  d.  Cells  were  then 
incubated  with  10  |jlM  BrdU  for  2  h.  Cell  were  fixed  with  1  ml 
of  4%  paraformaldehyde  in  PBS  and  then  permeabilized  with 
0.5%  Triton  X-100  in  PBS  for  5  min.  Cells  were  then  washed 
with  PBS  and  incubated  in  2  N  HC1  at  37° C  for  15  min  to 
denature  chromatin.  After  washing  with  PBS,  cells  were  incu¬ 
bated  with  2%  goat  serum  at  4°  C  for  18  h  to  block  nonspecific 
binding.  Cells  were  then  labeled  by  incubating  with  mouse  BrdU 
monoclonal  antibody  (Cell  Signaling)  for  2  h,  followed  by  incu¬ 
bating  with  goat  anti-mouse  Alexa  555  (Molecular  Probes). 
DAPI  (Life  Technologies)  was  used  to  stain  cell  nucleus. 

Colony  formation  assay 

Cells  were  transfected  with  miR-449a  mimic  or  control  oligos 
(25  nM)  for  24  h.  1000  cells  were  then  re-seeded  on  each  10  cm 
dish.  After  14  d,  colonies  were  visualized  by  staining  with  1% 
crystal  violet.  Colonies  were  counted  using  ImageJ  (NIH),  and 
differences  were  assessed  by  2-tailed  t-test. 


Construction  of  the  expression  vectors 

The  protein-coding  regions  of  the  MLAP4,  PKP4,  TSEN 1 5  and 
LEL1  genes  (without  the  3’UTRs,  to  avoid  the  down-regulation  by 
miR-449a)  were  amplified  by  PCR  from  human  cDNA  clones  (GE 
Dharmacon)  of  the  corresponding  genes.  The  amplified  MFAP4, 
TSEN  15  and  LEF1  sequences  were  inserted  into  the  EcoRI  and 
Xhol  restriction  sites  of  the  multiple  cloning  site  of  expression  vec¬ 
tor  pcDNA3. 1  +.  The  amplified  PKP4  was  inserted  into  the  Hindlll 
and  Notl  restriction  sites  of  pcDNA3. 1 +.  All  the  inserted  sequences 
were  verified  by  sequencing. 

miRNA  target  prediction 

miR-449a  target  sites  in  3rUTRs  of  the  mRNAs  were  identi¬ 
fied  based  on  seed  sequence  complementarity  using  IPA,  which 
identifies  any  7-nucleotide  region  (3'-  5')  in  a  given  3rUTR 
completely  complementary  to  the  seed  sequence  of  a  miRNA 
(nucleotides  2nd-8th  or  1st  -  7th,  5'-  3')  as  a  potential  target  site. 

mRNA  and  miRNA  expression 

Total  RNA  was  isolated  as  previously  described.74  mRNA 
expression  profiling  was  performed  using  the  Illumina  mRNA 
WG-6  v3  microarray  platform.  qPCR  measurement  of  mRNA 
expression  was  done  using  Taqman  Gene  Expression  Assay  (Life 
Technologies)  with  18s  rRNA  expression  as  a  control  for  normal¬ 
izing  RNA  loading.  miRNA  expression  was  measured  by  qPCR 
using  Taqman  micro  RNA  Assays  (Life  Technologies)  with 
expression  of  RNU44  RNA  used  as  a  loading  control. 

Luciferase  reporter  assay 

The  luciferease  reported  assay  was  performed  as  previously 
described.42  Briefly,  the  segments  of  the  wildtype  3/UTRs  con¬ 
taining  the  predicted  target  sites  of  miR-449a  were  cloned  from 
human  genomic  DNA.  Mutant  constructs  were  generated  with 
the  seed  target  sequences  mutated  to  CAGACGGA/U.  The 
3rUTRs  were  cloned  downstream  of  the  firefly  luciferase  coding 
sequences  into  the  pmirGLO  dual-luciferase  reporter  (Promega), 
a  vector  containing  both  firefly  and  Renilla  luciferase  cDNAs 
under  the  control  of  separate  promoter/terminator  systems.  The 
firefly  luciferase  was  used  as  the  primary  reporter  for  miRNA  reg¬ 
ulation  of  the  3’UTR.  The  Renilla  luciferase  is  an  internal  con¬ 
trol  for  normalization.  BE(2)-C  cells  were  co-transfected  with 
luciferase  reporters  (0.8  ng/ul)  and  miR-449a  mimic  or  control 
oligo  (25  nM).  Luciferase  activities  were  measured  after  3  d  using 
the  Dual-Glo  Luciferase  Assay  System  (Promega).  Firefly  lucifer¬ 
ase  activity  was  normalized  to  Renilla  luciferase  activity  to  evalu¬ 
ate  the  effect  of  the  miRNAs. 

Biotinylated-miR-449a  pulldown  assay 

The  pulldown  was  performed  using  a  protocol  based  on  Lai  A, 
et  al.46  The  protocol  and  the  calculation  for  the  target  enrich¬ 
ment  by  Bio-miR-449a  pulldown  is  illustrated  in  Suppl.  Fig¬ 
ure  7.  Briefly,  cells  were  transfected  with  either  Bio-miR-449a  or 
control  oligo  (Bio-Control)  (25  nM).  After  24  h,  cells  were 
lysed,  and  cell  lysate  was  separated  into  2  parts:  95%  was  subject 
to  biotin  pull-down  with  streptavidin-conjugated  beads,  and 
RNA  was  then  isolated  using  the  RNeasy  Kit  (Qiagen);  5%  was 
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used  to  isolate  RNA  as  an  input  control.  The  RNA  samples  were 
analyzed  by  qPCR  for  the  levels  of  individual  target  mRNAs 
with  GAPDH  mRNA,  which  is  not  a  target  of  miR-449a,  as  an 
internal  control  for  normalizing  the  target  mRNA  levels  in  each 
sample.  The  mRNA  enrichment  by  Bio-miR-449a  or  Bio- 
Control  pulldown  was  then  determined  by  the  ratio  of  mRNA 
levels  associated  with  the  pulldown  RNA  versus  the  input  RNA. 

Neuroblastoma  patient  survival  analysis 

Survival  analyses  were  based  on  2  published  neuroblastoma 
patient  cohorts:  the  Versteeg  cohort  which  includes  88  patients, 
and  the  Kocak  cohort  which  includes  476  patients.47,48  The 
microarray  gene  expression  and  the  patient  survival  data  for 
the  Versteeg  cohort  were  directly  obtained  from  http://r2.amc.nl. 
The  microarray  gene  expression  data  for  the  Kocak  cohort  were 
downloaded  from  http://r2.amc.nl,  and  patient  prognosis  data  for 
the  Kocak  cohort  were  downloaded  from  the  authors’  previous 
publication.48  For  patient  survival  studies,  overall  survival  was 
defined  as  the  time  from  diagnosis  until  death  or  until  last  contact. 
For  the  analysis  of  genes  in  each  of  the  patient  cohort,  the  patients 
were  divided  into  2  groups  based  on  whether  the  tumor  gene 
expression  level  of  the  patient  is  lower  (low  group)  or  higher  (high 
group)  than  the  average  gene  expression  of  the  patient  cohort.  Sur¬ 
vival  analyses  were  performed  using  the  Kaplan  and  Meier 
method,  and  significance  of  the  survival  difference  between  the  2 
groups  was  determined  using  2-sided  log-rank  test. 

Statistical  analysis 

The  statistical  significance  of  the  difference  in  frequency  dis¬ 
tribution  between  the  expression  of  the  predicted  miR-449a 


targets  and  non-targets  were  determined  by  2-sample  Kolmo- 
gorov-Smirnov  test,  with  P  <  0.05  considered  as  significant.  For 
evaluating  effect  of  the  panel  of  25  siRNAs  on  neurite  outgrowth, 
the  ^-values  for  neurite  lengths  associated  with  siRNAs  was  deter¬ 
mined  by  comparison  with  neurite  length  associated  with  the 
negative  control  oligo  using  multiple-sample  t-test.  We  consider 
a  siRNA  with  False  Discovery  Rate  (FDR)  <0.01  as  significantly 
inducing  neurite  outgrowth.  For  all  other  conditions,  the  statisti¬ 
cal  significance  for  each  treatment  was  determined  by  t-test  by 
comparing  the  treatment  group  with  control,  with  P  <  0.05  con¬ 
sidered  statistically  significant. 
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